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We present a study of the radiative decays of X(3872) at Belle. In 
the Xci7 final state, we got the first evidence of a new particle at 3823 
MeV/c 2 . 
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1 Introduction 



The Belle detector is a general purpose spectrometer built to test the Standard Model 
mechanism for CP- violation in B decays to charmonium (golden channel) pQ. Parallel 
to this, Belle has proven to be an ideal place to carry out charmonium spectroscopy 
thanks to its very clean environment. Many new cc and cc-like exotic candidates 
states such as rj c (2S), X(3872), X(3915), Z(3930), X(3940), Zi(4050)+ Z 2 (4250)+ 
F(4260), Z(4430) + and y(4660) have been found. Among them X(3872) is the most 
interesting state which was first observed in B^ — > (J/ipTT + iT~)K ± at Belle [2J. Soon 
after its discovery, it was confirmed by the CDF [3], DO [I] and BaBar j5] collabora- 
tions. Recently, it has also been observed at LHCb [6 J and CMS [7] . The observation 
of X(3872) in the same final states (J/if)7i + n~) in six different experiments, reflects 
the eminent status of the X(3872). X(3872)'s narrow width and the proximity of its 
mass, 3871.7 ± 0.2 MeV/c 2 to the D*°D° threshold makes it a good candidate for a 
DD molecule [8]. Other possibilities have also been proposed for the X(3872) state, 
such as tetraquark [9], ccg hybrid meson [10] and vector glueball models [TT] . 

Radiative decays of X(3872) provide a unique opportunity to understand the 
nature of X(3872). For example, X(3872) — > J/ip'y resulted in the confirmation of C- 
even (C = +) parity for X(3872) [T2J \T5\ Q3] . A similar decay mode X(3872) ->■ ^7 
can help to identify X(3872) as a charmonium, molecular or mixture of those states 
[HI [151 HE]- Belle reported no significant signal [14J and it disagrees with the evidence 
at BaBar p2]. 

A recent search for the charged tetraquark partner of X(3872)(— ^ J/^7r + 7r°) gave 
negative results [IT]. But still it is hard to totally rule out X(3872) as a tetraquark, as 
some tetraquark models predict X(3872) + to be broad, thus still difficult to observe 
at current statistics [18]. On the other hand, in both molecular and tetraquark 
hypothesis a C-odd parity (C = — ) partner can exist and it may dominantly decay 
into Xcil and \c2l final states. Partner searches would be another approach to reveal 
X(3872) structure. Along with this, undiscovered 3 D 2 charmonium (f/> 2 ) is expected 
to have a significant branching fraction to Xcil and Xc2l [T9~| BP] . In this report, we 
also describe the first evidence of ip2 — >■ Xal decay. 

2 Reconstruction 

B mesons reconstructed using B ± — > (Xci( — > J l' l Pl)l)K ± &n.&B ± — > (Xc2(— > J '/4'l)l)K ± 
decay modes are used in the search for C-odd partner of the X(3872) and other new 
narrow resonances. The results presented here are obtained from a data sample of 
772 x 10 6 BB events collected by the Belle detector pi] at the KEKB [22] energy- 
asymmetric e + e~ collider operating at the T(4S') resonance. 

The J/ip meson is reconstructed in its decays to £ + £~ (£ = e or fi). In e + e~ decays, 



1 



the four-momenta of all photons within 50 mrad of each of the original e + or e~ tracks 
are included in the invariant mass calculation [hereafter denoted as M e + e -( 7 \], in order 
to reduce the radiative tail. The reconstructed invariant mass of the J/ip candidates 
is required to satisfy 2.95 < M e+e - (7) < 3.13 GeV/c 2 or 3.03 < M^+^ < 3.13 GeV/c 2 . 
A mass- and vertex-constrained fit is applied to all the selected J/ip candidates in 
order to improve their momentum resolution. The Xci,c2 candidates are reconstructed 
by combining J/ip candidates with a photon having energy (E^) greater than 200 
MeV. The photons are reconstructed from the energy deposition in electromagnetic 
calorimeter. To reduce the background from 7r° — > 77, we calculate a likelihood 
function to distinguish an isolated photon from 7r° decays using the photon pair's 
invariant mass, smaller photon's laboratory energy and polar angle [23]. Then we 
reject a photon having the 7r° likelihood ratio greater than 0.7 by combining with any 
other photon. The reconstructed invariant mass of xa and Xc2 is required to satisfy 
3.467 < M J/i>y < 3.535 GeV/c 2 and 3.535 < M Jf ^ < 3.611 GeV/c 2 . A mass- and 
vertex-constrained fit is again performed to all the selected Xci and Xci candidates in 
order to improve their momentum resolution. Charged tracks are identified as a kaon 
using information from the particle identification devices. 

To reconstruct the B candidates, each Xci, Xc2 is combined with a kaon candidate 
and an additional photon having E^ > 100 MeV. In this additional photon selection, 
we reject the photon which in combination with another photon in that event, gives 
mass in the region around ir° mass defined as 117 < M 77 < 153 MeV/c 2 . In order to 
remove the reflection of photons coming from Xci and Xc2 decays, we reject the best 
Xd or Xc2 daughter photon from the additional photon list to form a B candidate. 
To identify the B candidate, two kinematic variables are used : energy difference 
AE = Eg — El eam and beam-energy constrained mass M bc = J (E£ eam ) 2 — (p^" s ) 2 , 
where E* B is the center-of-mass frame (cms) beam energy, and E* B and p* B are the cms 
energy and momentum of the reconstructed particles. In case of multiple candidates, 
AE closest to is selected as the best one. Invariant mass of the final state (M Xcl7 
and M Xc27 ) is used to identify the resonance. In order to improve the resolution of 
M XciCI jMj, we scale the energy of 7 to make AE equal to zero. 

To suppress continuum background, events having a ratio of the second to zeroth 
Fox- Wolfram moments [25J R2 > 0.5 are rejected. Large B — > J/ipX MC samples 
(corresponding to 100 times the data sample size used in this analysis) are used to 
study the background. The non- J/tp (non-x CX7 ) background is studied using the 
M u (Mj/^) sidebands in data. B ± — > (x&l)K ± and B ± — > (xc2l)K ± yields are 
extracted from a 2D UML fit applied to the distribution in the M Xcx7 - Mt, c space. 
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Figure 1: 2D UML fit projections in M Xcl7 (top) and M Xc21 (bottom) for the signal 
region M^ c > 5.27 GeV/c 2 . The curves show the signal [red large-dashed for ip', pink 
dashed for ip 2 and yellow dotted for X(3872)] and the background component [purple 
dotted-dashed for combinatorial, dark green two dotted-dashed for B ± — > -0' (other 
than Xcil)^^ and cyan three dotted-dashed for peaking component] as well as the 
overall fit [blue solid]. 

3 Results 

In B ± — > {xcil)K ± decay, a clear peak of ip' — > Xcil is observed in the Xcil invariant 
mass (M Xcl7 ) projection. In addition to this, we also find clear evidence of a narrow 
peak at 3823 MeV/c 2 , denoted as X(3820) hereafter. Since statistics are still limited, 
its width is poorly constrained to be 4 ±6 MeV even if it is floated in the fit. Therefore 
we set X(3820) natural width to be MeV for further discussion. We estimate the 
statistical significance of X(3820) — > Xcil to be 4.2<r including systematic uncertainty. 
In our search for X(3872) — > Xcil, no signal is seen and 90% confidence level (C.L.) 
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Channel 


Yield 


£(10- 4 ) 


B ± -> XciT)^ 
B± _> ^ Xc27 )K± 


193.21^ 


7 7+u.a+o.y 

' • ' -0.7-0.8 

6.3 ±0.9 ±0.6 






xl0~ 6 


5± ->X(3820)(-> XciT)^ 
5 ± ^X(3872)(^ XciT)^ 
S ± ^X(3820)(^ Xc2l)K ± 
£±^X(3872)(^ XcaT)^ 


33.21^ 
-0.9 ±5.1 

0.311:? 

4-71^ 


q 7 +a.8+i.i 

J - ' -2.5-1.0 

< 2.0 

< 3.6 

< 6.7 



Table 1: Summary of the results. Measured B (with 90% confidence level (C.L.) 
upper limit (U.L.) for B ± ->■ X(3872)(^ Xcff)^, B± ~* ^(3823) (-> Xc2l)K ± and 
B ± ->■ X(3872)(^ Xc2l)K ± decay modes). B for £ ->■ XX is B(S -> XK)B{X -> 
X cx 7), here X stands for X(3820) and X(3872). For 23, the first (second) error is 
statistical (systematic). 



upper limit (U.L.) on B{B ± -> X(3872) J FT ± )£(X(3872) -> Xci7) is estimated using a 
frequentist approach. In B ± — > (xc2j)K ± , we observe a clear peak of the ip' — > Xc2l- 
However, we do not see any hint of a narrow resonance with the current statistics. 
Figure [TJ shows the projections of 2D UML fits to the M Xcl7 and M Xc27 cases for the 
signal region, M bc > 5.27 GeV/c 2 . Tabled] summarizes the results of B ± — » {xcil)K ± 
and 5 ± — > (xc2l)K ± decays. 

The mass of X(3820) is near the potential model expectations for the center-of- 
gravity (cog) of l 3 Dj states: the Cornell [26] and the Buchmuller-Tye [27] potential, 
which gives M cog (ID) = 3810 MeV/c 2 . Some models predict the mass of 3 D 2 (J PC = 
2~ ) states to be 3815-3840 MeV/c 2 HUES]. X(3820) mass agrees quite well within 
the expectation. The X(3820) state is likely to be the missing 3 D 2 cc (ip 2 ) state. 

4 Summary 

In the study of B ± — >■ (xdl)K ± and B ± — > (xc2l)K ± decays, Belle observe the ip' 
peak in both decay modes and other peak is found to be consistent with expectation. 
In B ± — > (xcil)K ± 1 Belle find the first evidence of a narrow state having mass of 
3823.5 ± 2.5 MeV/c 2 . This narrow state is likely to be the missing ip 2 { 3 D 2 cc state) 
because the observed mass totally agrees with a theoretical expectation and Xci7 is 
one of the dominant decay modes. While in the search of C-odd partner of X(3872), 
no signal is found at the current statistics and U.L. in 90% C.L. are provided. 
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